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Abstract 

The detection of acoustic signals from ultra-high energy neutrino interactions is a promising method to measure 
the flux of cosmogenic neutrinos expected on Earth. The energy threshold for this process depends strongly on the 
absolute noise level in the target material. The South Pole Acoustic Test Setup (SPATS), deployed in the upper part of 
four boreholes of the IceCube Neutrino Observatory, has monitored the noise in Antarctic ice at the geographic South 
Pole for more than two years down to 500 m depth. The noise is very stable and Gaussian distributed. Lacking an 
in-situ calibration up to now, laboratory measurements have been used to estimate the absolute noise level in the 10 
to 50 kHz frequency range to be smaller than 20 mPa. Using a threshold trigger, sensors of the South Pole Acoustic 
Test Setup registered acoustic events in the IceCube detector volume and its vicinity. Acoustic signals from refreezing 
IceCube holes and from anthropogenic sources have been used to test the localization of acoustic events. An upper 
limit on the neutrino flux at energies E v > 10 1 1 GeV is derived from acoustic data taken over eight months. 

Keywords: acoustic neutrino detection, absolute noise level, neutrino flux limit 
PACS: 43.58,+z, 43.60.+d, 93.30.Ca 



1. Introduction 

During recent years it has been extensively studied, 
whether the glacier ice at the South Pole is a suitable 
material for the detection of cosmic neutrinos above 
10 I8 eV energy, using their acoustic signals emitted in 
the 10 to 50 kHz region. With data from the South Pole 
Acoustic Test Setup [1], the speed of sound has been 
measured down to 500 m depth. It was found to be con- 



stant below 200 m, favorable for neutrino detection Q . 
The sound attenuation in the 200 to 500 m depth region 
was measured to be much stronger than expected fl. 
An important precondition for acoustic neutrino detec- 
tion is a precise knowledge of the steady and transient 
noise levels on top of which a possible neutrino signal 
has to be detected. 

Long duration noise studies in deep water of the 
Mediterranean Sea and Lake Baikal JfJ have 
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shown strong variations with time due to changing en- 
vironmental conditions such as wind speed and rain. 
Sound sources such as ships and animals (e.g. sperm 
whales) contribute as well. In quiet periods, however, 
noise levels as low as a few mPa are reported for the 
frequency region of interest. 

Acoustic noise at the South Pole is assumed to be low 
and stable because none of the sources mentioned above 
are expected to contribute. Due to the depth dependence 
of the density, the first 200 m of firn ice acts as an acous- 
tic filter. Anthropogenic and environmental noise enter- 
ing the ice from the surface will be refracted back. 

The glacier moves with a velocity of about 10 m per 
year over the bedrock beneath. Whether this or other 
phenomena give rise to ultra-sound of detectable ampli- 
tude is one subject of the investigations presented in this 
paper. 

In Sec. [2] the detector geometry, calibration and the 
various data taking modes will be described. In Sec. [3] 
the permanent noise conditions will be discussed, fol- 
lowed by information about the detection of sound sig- 
nals produced in the ice by construction activities of the 
IceCube experiment [7] in Sec.|4] Finally, a limit for the 
flux of cosmic neutrinos will be derived from acoustic 
data accumulated during eight months. 

2. Detector setup and data taking 

2.1. Detector geometry 

To measure the acoustic properties of the ice at the 
geographic South Pole, the South Pole Acoustic Test 
Setup (SPATS) Ht], a system of four instrumented verti- 
cal lines, called strings A, B, C, and D, was installed 
in boreholes of the IceCube neutrino observatory af- 
ter deployment of the IceCube optical modules. Each 
string holds seven stages, referred to as stages 1 to 7 
from top to bottom; each stage is a combination of an 
acoustic sensor and a transmitter, vertically separated 
by ~ 1 m. The horizontal distances between strings 
range from 125 m (which is the spacing between Ice- 
Cube strings) to 543 m. Vertically, the depth range from 
80 to 500 m is instrumented with increasing spacing of 
sensors in the deeper ice to be able to sample the transi- 
tion from the firry region to the bulk ice. On strings A, 
B, and C the stages are located at depths of 80, 100, 140, 
190, 250, 320, and 400 m; string D is deployed deeper 
in the ice and the stages are located at 140, 190, 250, 
320, 400, 430, and 500 m depth. 



3 The transition region from a snow/air mixture at the surface to 
solid ice is called firn. It has a width of about 170 m at South Pole Q|. 



The SPATS sensors consist of a cylindrical stain- 
less steel pressure housing with a diameter of 10 cm 
in which three piezoelectric elements are mounted to 
the wall, separated by 120 degrees to provide full az- 
imuthal coverage. A three-stage low noise pre-amplifier 
with a gain of 10 4 is attached directly to each piezo- 
electric element and the analog signal is transmitted to 
the surface via a shielded twisted pair cable. On string 
D at depths of 190 m and 430 m, an alternative type of 
sensor, HADES (Hydrophone for Acoustic Detection at 
South Pole), was deployed. For the HADES sensors, the 
piezoelectric element and pre-amplifier have been cast 
in resin and mounted below the housing. This allows 
us to study acoustic signals with different systematics 
introduced by the sensor. 

The sensor channels are each identified by a label of 
the form XSy(z), where X is the letter of the string, S 
refers to a sensor, y is the number of the stage, and z the 
channel number within the sensor (running from to 2). 

2.2. Data acquisition 

On the surface, a read-out box buried in the snow 
above each string contains an industrial PC, called 
string-PC, that is used for digitization, time stamping, 
and storage of the data. Each string-PC is connected by 
a symmetric DSL connection to the SPATS master-PC 
that is housed in the IceCube laboratory at the South 
Pole station. The master-PC collects the data from all 
four string-PCs, distributes a GPS timing signal to them, 
and prepares the data for transfer to the northern hemi- 
sphere via satellite or for tape storage. 

Since 28 August 2008, SPATS has been operating as 
a detector for transient acoustic signals. Out of each 
hour, 45 minutes are used for triggered data taking while 
in the remaining 15 minutes environmental monitoring 
and system health data are recorded, including measure- 
ments of the noise floor. Details on the DAQ system can 
be found in Ref . [ 1 ] . To monitor the noise floor, every 
hour the three channels of each sensor simultaneously 
record 100 ms of untriggered data with a sampling rate 
of 200 kHz. In transient data taking mode, three sensor 
channels from each string are used, located at three dif- 
ferent depths. The data stream from each channel, dig- 
itized with a sampling rate of 200 kHz, is continuously 
monitored in the string-PC. A sensor noise band is de- 
fined for each channel i spanning symmetrically around 
the channels noise mean value from /y, - 5.2 cr, to 
ju, + 5.2 <x, as illustrated in Fig. [2] where cr, is the width 
of the noise distribution of the respective channel. If an 
ADC value x, outside the noise band (x, > ju, +5.2 cr, or 
Xi < jU, -5.2 cr, is measured, 5 ms of data (1001 samples) 
are written to disk, centered on the triggering sample. 
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The threshold of 5.2 cr, was chosen to fit the amount 
of data into the allocated daily satellite bandwidth of 
150 MB. The detector has been operated in two different 
modes, hereafter called mode 1 and mode 2. In mode 1, 
the sensors located at 190, 250, and 320 m depth were 
used for transient data taking. In mode 2 a deeper sensor 
configuration (250, 320, 400 m) was chosen. 

One transient data file per string is written every hour. 
Transient data taking runs are successful in 93% of all 
cases. The remaining runs commonly fail due to DAQ 
problems. No correlation between failing runs and large 
amplitude pulses, which would bias the analysis de- 
scribed in this work, has been observed. During mea- 
surements with the acoustic "pinger" J2, 0] large am- 
plitude pulses have been recorded with a rate of up to 
10 Hz without an increase in the occurrence of failed 
runs. 

As can be seen in Fig. |l(a)| there are quiet periods 
during the austral winter and more noisy periods during 
summer with IceCube hole drilling and string deploy- 
ment taking place. We distinguish between four periods 
of data taking; details can be found in Tab. \T\ On av- 
erage, about 65000 triggers are collected per day. The 
data are dominated by sensor noise, i.e. triggers caused 
by single ADC samples from the tails of the Gaussian 
noise distribution. In the quiet periods, these account 
for 99.99% of all recorded waveforms, for which no co- 
incident signals in other sensor channels are observed. 
Nearly all of these triggers do not contain typical wave- 
forms but have only a single spike larger than the trigger 
threshold. 

2.3. Sensor calibration 

The sensitivitjQ and equivalent self-nois^f] of all 
SPATS sensors were measured in the laboratory prior 
to deployment using the comparison method. The mean 
sensitivity (including the pre-amplifier) averaged over 
all sensor channels and all frequencies in the relevant 
frequency band from 10 to 50kHz is 2.8 + 0.8 V Pa" 1 
(equal to -llldBre. 1 VyuPa 1 ), where the error in- 
dicates the spread between the different sensors. The 
equivalent self-noise in the same frequency band, aver- 
aged over all SPATS sensors, is 7 mPa. Details of the 
measurement are presented in Ref. [1]. The calibration 
was performed in water at a temperature of 0°C and at 



Sensitivity is defined as the ratio V/p of the voltage V induced at 
the output of an acoustic sensor, given a pressure amplitude p, with 
units of VPa -1 . 

5 Equivalent self-noise is defined as the electronic self-noise of the 
device converted into an equivalent sound pressure level using the sen- 
sor's sensitivity. 



normal pressure. It is not clear a priori, whether the 
sensitivity of the sensors remains unchanged when de- 
ployed deep in Antarctic ice at a temperature of -50°C 
and an increased but not well known static ambient 
pressure. The water in the IceCube boreholes starts to 
freeze from the top and from the sides so that cavities 
are formed in which the static pressure will increase 
above the pressure of the water column at sensor depth. 
After the hole is completely frozen, relaxation of the 
freshly frozen "hole ice" to the surrounding bulk ice oc- 
curs on an unknown time scale. Furthermore, the sen- 
sitivity of the sensor can be influenced by the different 
acoustic coupling, determined by the different acoustic 
impedance matching between the water-sensor interface 
and the ice-sensor interface, by the unknown horizontal 
position of the sensor within the IceCube hole, and by 
possible shadowing of parts of the sensor by the Ice- 
Cube main cable running down the hole parallel to the 
sensors. 

Due to the absence of standardized acoustic sources 
for ice, it is not possible to recalibrate the sensors after 
deployment. Also, it is beyond the scope of the SPATS 
project to reproduce the combined influence of very low 
temperatures, high static pressure and ice in a labora- 
tory environment. However, we can study the influence 
of these aspects on the sensitivity by separate investi- 
gations. In the following discussion, we will assume 
that the influences of the environmental parameters are 
uncorrected and can thus be used to estimate the in-situ 
sensor sensitivity from the results of different laboratory 
experiments. A typical SPATS sensor was calibrated in 
air at different temperatures. A linear increase of its sen- 
sitivity with decreasing temperature was measured. The 
sensitivity increases by a factor of 1 .5 + 0.2 from 0°C 
to -50°C yj]. The same sensor was calibrated at room 
temperature in a water filled pressure vessel at differ- 
ent static pressures from 1 bar to 100 bar. No systematic 
change in the sensitivity with pressure was observed; 
the sensitivity was found to be stable within 30% 01 • 
The effect of the acoustic coupling of the sensor to the 
ice, which can lead to a flattening of the frequency re- 
sponse function and a reduction of the sensitivity due to 
a damping of mechanical resonances of the sensor hous- 
ing, is not yet understood and will be studied in an ice 
block in the laboratory using the reciprocity calibration 
method (for a discussion of the reciprocity calibration 
method, see e.g. Ji)]). 

Since the correlation between the different environ- 
mental parameters is unknown when the sensor is ex- 
posed to a combination of all of them in the deep 
Antarctic ice, the determination of the absolute noise 
level is currently only roughly possible. Multiplying the 
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Figure 1; (a) Number of acoustic hits per day and (b) acoustic events with hits at all 4 strings since 1 January 2009. Drill periods are indicated in 
light gray. For the definition of a "hit" and an "event" see Sec. [4] 



Table 1 : Characteristics for different data taking periods of transient noise triggers. 



Name 


Quiet period 1 


Drill period 1 


Quiet period 2 


Drill period 2 


Stall date 


28 Aug. 2008 


1 Nov. 2008 


1 Mar. 2009 


1 Nov. 2009 


Duration (days) 


65 


120 


245 


120 


Available files 


5820 


9845 


22664 


10567 


Avail./total 


0.93 


0.85 


0.96 


0.92 


Detector mode 


1 


1 


2 


2 
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sensor's sensitivity changes due to temperature (xl.5) 
and due to pressure (xl) and adding the uncertainties in 
quadrature, we find that the sensitivity of the sensor in 
the deep ice is increased by a factor of 1 .5 + 0.4 as com- 
pared to the value obtained in the laboratory, resulting 
in a mean sensitivity of 4.2 + 1.6 V Pa -1 . 

3. Properties of the noise floor 

3.1. General 

We observe a Gaussian distribution of ADC values 
for each sensor channel. Thus the noise can be char- 
acterized by two parameters: a mean value and a stan- 
dard deviation. The mean value depends on an instru- 
mental DC offset in the different channels and is always 
close to zero. The standard deviation is a measure for 
the noise level in the sensor, which is a superposition 
of sensor electronic self-noise, electromagnetic interfer- 
ence picked up on the signal cable from the sensor to the 
surfac^, and possible acoustic noise contributions from 
the surrounding ice. 

The distribution of ADC values from 73.1s of 
recorded noise data is shown in Fig. [2] The dynamic 
range of the 12-bit ADC is ±5V, corresponding to 
2.4 mV per ADC count. The data are perfectly de- 
scribed by a Gaussian. The mean value and trigger 
thresholds for transient data taking are indicated in the 
graph. The four samples outside the noise band match 
very well the expectation from the average SPATS trig- 
ger rate of 4.6 triggers during these 73.1 s. 

3.2. Stability 

We have monitored the noise level in all sensor chan- 
nels for more than three years, beginning with the de- 
ployment of the first SPATS sensors in January 2007. 
Figure [3] shows the RMS of the noise as a function of 
time for three typical sensor channels on string C that 
participate in the transient data taking. All available 
data from deployment until autumn 2010 is shown. It 
can be seen that, apart from some short-time excesses 
that will be discussed below, the noise level is very sta- 
ble, the typical fluctuations being cr RMS / (RMS) < 10~ 2 . 

In 2007, the first year of SPATS operation, we mea- 
sured a higher and less stable noise level. During this 
period, the sensors were powered on only for data tak- 
ing and powered off afterwards which causes them to 



Electromagnetic interference is expected to be small since the sig- 
nal is transmitted differentially from the sensor to the ADC on the 
surface. 
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ADC value 

Figure 2: Distribution of ADC values from all noise data recorded for 
channel CS7(2) in July 2009 (73.1 s of data in total). The red dashed 
line is a Gaussian fit to the data; the vertical lines indicate the mean 
value p (dashed) and the trigger thresholds at )i ± 5.2 cr. The peak at 
ADC value is an understood feature of the ADC. 
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Figure 3: Noise RMS (from zero to 100 kHz) as a function of time for 
three sensors on string C participating in transient data taking. The 
time window spans all data from deployment till today. 
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Figure 4: Noise RMS as a function of time during freeze in for sensor 
channel DS7(2) deployed at 500 m depth. 



heat up during the measurement and change their self- 
noise characteristics. Since December 2007, all sensors 
are powered on continuously and are in thermal equilib- 
rium with the surrounding ice. Noise studies following 
a power outage show that it takes several hours for a 
sensor to reach thermal equilibrium after it is powered 
on. The short-term noise excesses, which occur only in 
the Austral summer seasons, can be correlated to Ice- 
Cube deep-ice drilling activity. The visible spikes cor- 
respond to the holes drilled closest to the SPATS array. 
Due to technical reasons, data on the noise level dur- 
ing the freeze in of the sensors after deployment is only 
available for string D which was installed one year af- 
ter the other three strings. It is shown for one sensor 
channel in Fig. [4] 

We observe an increase of the noise level after the de- 
ployment of the sensor on 24 December 2007 that lasts 
for about three weeks, after which the noise level be- 
came stable. On top of that, excesses correlated with 
IceCube deep-ice drilling can be seen. We interpret the 
rise of the noise level as a combination of the increase 
of sensitivity with decreasing temperature (cf. Sec. [23} 
and an improved acoustic coupling of the sensor to the 
bulk ice. 

3.3. Determination of the absolute noise level 

To determine the absolute noise level from the data, 
i.e. the sound pressure incident on the sensor, the sen- 
sitivity needs to be known. In general, the sensitivity 
will be a function of direction and frequency. In the 
absence of a specific noise source model, an equivalent 
acoustic power at the position of the piezoelectric ele- 
ment is derived assuming isotropic noise. The power is 
then translated into an effective pressure amplitude. 



Channel CS5(1) (250m depth) 




90 100 
Frequency (kHz) 



Figure 5: Average power spectral density (PSD, shown as solid lines) 
and distribution for each frequency bin (gray scales) for three sen- 
sors on string C participating in transient data taking. All noise data 
recorded in July 2010 are shown. For comparison, the sensitivity of 
the sensor channels measured in the laboratory prior to deployment is 
shown as dashed lines. 



In Fig. [5] we show the voltage power spectral density 
(PSD) distribution for three sensor channels participat- 
ing in transient data taking. The plot was obtained as 
follows: all noise floor data recorded in July 2010 were 
Fourier transformed in sets of 1000 samples each (fre- 
quency resolution A/ = 200 Hz) and the correspond- 
ing PSD values were filled into a two dimensional his- 
togram. The gray scale is a measure for the probabil- 
ity of occurrence of a certain PSD value at a given fre- 
quency. The solid lines represent the mean values, cal- 
culated on a linear PSD scale, in each frequency bin. 
The error of the mean is also indicated, but too small to 
be visible. 

Figure [5] demonstrates that the spectral shapes differ 
between all sensors. The sensitivity of the sensors as a 
function of frequency, as measured in water in the labo- 
ratory prior to deployment, is indicated as dashed lines. 
It is expected that the resonance structure of the PSD is 
mainly governed by the mechanical response of the sen- 
sor and not determined by the spectrum of the incident 
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acoustic background noise on the sensor. Especially the 
steel housing and its coupling to the piezoelectric ce- 
ramic, which can be slightly different for different sen- 
sors, should have a large effect. It can be seen in Fig. [5] 
that peaks in the sensitivity are not reflected as peaks 
in the PSD as would be expected for a smooth acoustic 
noise spectrum. This supports the assumption that the 
resonance behavior of the sensor, and thus its sensitiv- 
ity, is modified by the coupling of the sensor housing 
to the ice. Due to the suspected change in the spectral 
sensitivity during freeze-in, we will not calculate the ab- 
solute noise level by dividing the PSD by the sensitiv- 
ity to determine the noise spectrum in units of pressure 
density and integrate over the relevant frequency range. 
This procedure would introduce unknown errors by un- 
derestimating the contribution from frequency regions 
with high sensitivity in the laboratory calibration. In- 
stead, we assume a single mean sensitivity for all sen- 
sors, determined by averaging the laboratory sensitiv- 
ity of each sensor over the frequency range from 10 
to 50 kHz and subsequently averaging over all sensors 
and applying the correction factor for temperature and 
pressure. This procedure yields a mean sensitivity of 
(S 10-50) = 4.2 ± 1.6 VPa~' as discussed in Section [2731 
We determine the absolute noise level from each sensors 
voltage PSD integrated from 10 to 50 kHz, using the 
July 2010 data presented in Fig. [5] Attenuation losses 
in the cable of -0.6 dB/ 100 m are corrected for. This 
assumes the worst case scenario that all the measured 
noise is produced in the sensor or is acoustic noise in the 
ice and no additional electromagnetic noise is induced 
during transmission. Noise induced further upstream in 
the DAQ chain would be over-corrected for cable at- 
tenuation and result in an overestimation of the noise. 
Figure [6] shows the resulting noise level for all opera- 
tive SPATS channels. We separate the sensors into two 
groups: sensors above 200 m depth and sensors below 
200 m. The latter ones are used for the transient noise 
analysis in the remainder of this work. For the shallow 
sensors, we calculate a mean noise level of 21 mPa with 
a 5 mPa (lcr) spread between the data points. The aver- 
age noise level in the deep sensors is (16 + 3) mPa. This 
still includes the contribution from electronic self-noise, 
that has been measured in the laboratory prior to de- 
ployment to be 7 mPa equivalent on average. Subtract- 
ing this contribution quadratically leads to an estimated 
mean noise level in South Polar ice of 20 mPa (shallow) 
and 14mPa (deep) integrated over the frequency range 
relevant for acoustic neutrino detection of 10 to 50 kHz. 
Using the simulation described in Sec.[5]and assuming 
an acoustic attenuation length of 300 m, a pressure of 
14 mPa corresponds to the amplitude of the acoustic sig- 
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Figure 6: Estimated absolute noise level integrated from 10 to 50kHz 
for all SPATS channels. The error bars indicate the uncertainty on the 
sensitivity of the sensor channels. Sensors above 200 m depth (stages 
1 to 4) and below 200 m (stages 5 to 7) are treated separately (see 
text for details). Sensors AS3 and CS1 are broken and no data are 
available. The solid lines indicate the mean values and the dashed 
lines indicate the lcr spread of the data points. An average equivalent 
self-noise of 7 mPa can be subtracted quadratically for all sensors. 



nal generated by a neutrino of energy 10 1 1 GeV interact- 
ing in a distance of about 1000 m. 

The origin and significance of the decrease of the 
noise level with depth, that is visible in Fig. [6] remains 
unclear. One possible qualitative explanation for the 
observed depth dependence is a contribution of noise 
generated on the surface. Due to the gradient in the 
sound speed with depth [2], all noise from the surface 
will be refracted back towards the surface, thus shield- 
ing deeper regions from surface noise. 



4. Transient noise events 

A triggered waveform ("hit") contains 500 samples 
measured before the trigger and 500 samples after the 
trigger, each separated by 5 jus (see also Sec. 12.2b . If the 
acoustic signal lasts longer than 2.5 ms, the following 
trigger is added to the hit under consideration. Hits from 
all strings are ordered in time offline and merged into 
one file per day. This file is processed through a clus- 
ter algorithm to find all consecutive hits within 200 ms, 
the time necessary for an acoustic signal to cross the 
SPATS array. All hits per cluster are considered to form 
an acoustic event. Events with more than 5 hits from at 
least 3 strings are used to localize their source position 
in the ice (see Section |4T1 ). The distribution in time of 
acoustic events with hits from all four strings is shown 
inFig.[T(b)l 
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4.1. Source vertex reconstruction 

The acoustic event reconstruction algorithm is based 
on the solution of the system of equations with n = 
1,...,4 

(x„ - xo) 2 + (y„ - y ) 2 + (z„ - zof = [v s (t n - t )] 2 . (1) 

Only four sensors with their signal arrival times and 
positions t n , x n ,y„,z n are used in a single reconstruc- 
tion. The calculated event vertex is located at the space 
time point to,xo,yo,zo, where the z-axis points verti- 
cally upward and z = corresponds to the ice sur- 
face. The signal velocity in ice is taken to be constant 
with v s = 3878 m/s J^L and the propagation direction 
is assumed to be straight. The assumption of a con- 
stant speed of sound is only suitable for events below a 
depth of around 200 m and leads to a spread of recon- 
structed event positions for shallower depths, as one can 
see from simulations (Section [43J. Solving the system 
of equations above provides an event vertex for a sin- 
gle four sensor combination. With twelve sensors in the 
used SPATS configuration, statistical predictions can be 
made by using all possible combinations i = 1, m of 
four sensors on four different strings per acoustic event. 
In case of a noise hit in a sensor, the reconstruction al- 
gorithm for this combination does not converge or the 
result lies far outside the sensitive SPATS area. The 
deviation from the mean vertex position of all possible 
sensor combinations 
1 m 

?=-2> (2) 

1=1 

with r t - (x' Q ,y' ,z' ), is used to improve upon the back- 
ground rejection. This is done by rejecting reconstruc- 
tion results for a single sensor combination if the dis- 
tance from the mean position in x,y or z is above 250 m. 

4.2. Acoustic event sources 

In Fig. |7(a)| all reconstructed four-string events are 
plotted according to their position in the IceCube co- 
ordinate frame. Figure |7(b)| shows their location with 
respect to the SPATS strings (large circles) and IceCube 
holes (small circles). The triangles reflect the positions 
at which a Rodriguez- Well (RW for short) [10] is lo- 
cated (see Table|2|. Such wells are used for the produc- 
tion and cycling of water for the IceCube hot water drill 
system. As can be seen, the acoustic events are concen- 
trated either at IceCube boreholes or at Rodriguez-Well 
locations. In Fig. |7(c)| the event depth distribution ver- 
sus time is shown. Almost no events are located above 
50 m depth. In quiet periods, events are concentrated 
between 80 and 150 m. During drilling periods vertices 
are still found down to 600 m. 



4.2.1. Events from IceCube holes 

Acoustic events were observed from nearly all Ice- 
Cube holes drilled in both seasons, when transient data 
taking was active. 

The statistics collected in the second season was 
much larger, due to the fact that most of the 2009/10 
holes were located close to the center of the SPATS- 
detector. The event distributions derived from different 
holes are very similar. We investigate 2093 events from 
hole 81, the hole with the highest statistic, as an ex- 
ample. Events are observed for 20 days in the hole re- 
gion (±20 m with respect to the center of hole 81) during 
the periods of firn ice drilling (< 50 m depth), bulk ice 
drilling (50 to 2500 m depth) and refreezing, as can be 
seen in Fig. |8(a)| Figure [8(b)| shows the depth distribu- 
tion of the events versus time. Before drilling, events are 
observed at 40 to 100 m depth probably connected with 
noise from the firn drill hole. During the procedure of 
hot water drilling a few related events are found. Strong 
sound production starts about three days after drilling 
is finished, due to the refreezing process. About 30% 
of the registered events from this hole are concentrated 
in two spots at 120 m and 250 m depth but reach down 
to about 600 m. The reason is that the hole does not 
refreeze homogeneously, but forms frozen ice plugs be- 
tween regions that are still filled with water. The pres- 
sure produced in this way may give rise to cracks near 
the ice water boundary which would appear with sound 
in the 10 to 100 kHz frequency region. Relaxation later 
continues within "arms" freezing towards the hole sur- 
face and down to the lower ice plug (see Fig. |8(b)| >. 
Besides providing information about the refreezing pro- 
cess of water filled IceCube holes, one can also use the 
corresponding acoustic events to understand the preci- 
sion of the vertex localization algorithm. In Fig. |9(a)| 
and Fig. |9(b)| the reconstructed x and y event positions 
are shown for hole 81 in the IceCube reference coor- 
dinate system. The average values including statistical 
uncertainties for the (x,y) position of hole 81 are de- 
termined to (42.0 + 0.1) m in x and (38.5 ± 0.1) m in 
y. The width of the distributions is 2.6 m and 5.0 m re- 
spectively, to be compared with a hole diameter of about 
0.7 m. The calculated values deviate from the actual 
hole positions at the surface by 0.4 m (in x) and 3.0 m 
(in yo The possible reason for this deviation will be 
discussed in the simulation section (Sec. [43} below. 



7 The actual positions of the sensors in the x-y plane in the holes are 
known with a precision of 0.5 m due to the hole width and inclinations 
versus depth from the drilling process. 
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Figure 7: Shown are (a) the relative abundance of reconstructed acoustic events in the horizontal plane of the IceCube coordinate system and (b) 
the actual vertex position of all transient events recorded since August 2008. The sources of transient noise are the Rodriguez Wells (RW), large 
caverns melted in the ice for water storage during IceCube drilling, and the refreezing IceCube holes. Dark gray circles (IceCube 2): positions of 
IceCube holes drilled in the period of transient data taking, light gray circles (IceCube 1): other IceCube holes, black filled circles: locations of 
SPATS strings, triangles and square: location of RW. The elongated structures are discussed in Sec. 14.31 In (c) the depth distribution of acoustic 
events versus time since August 2008 is shown. The zero at the time axis corresponds to Jan 1 , 2009 and drill periods are indicated in light gray. 
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Figure 8: (a) Number of acoustic events per hour from hole 81 before and after drilling and (b) depth distribution versus time of acoustic events 
from hole 8 1 . The vertical line indicates the end of hot water drilling for hole 8 1 . 



800r 



700: 



E 600: 



500 : 



400 : 



a 300- 
E 

i 200: 
100: 



350 



20 30 40 50 
x-coordinate in m 

(a) 



60 



70 




20 30 40 50 60 
y-coordinate in m 

(b) 



Figure 9: (a) x-coordinate and (b) y-coordinate distribution of acoustic events from hole 81. The vertical line shows the nominal position of hole 
81. 
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4.2.2. Noise from Rodriguez-Wells 

When first acoustic events had been reconstructed 
during the first period from August to November 2008 
(quiet period), a strong clustering in a certain region of 
the x-y plane at about (-150 m, 300 m) became visible. 
It was found that this was the position of the 2007/08 
Rodriguez- Well used for the hot water drilling system. 

This type of well has been introduced by Rodriguez 
and others in the early 1960s [10] to supply water from 
a glacier in Greenland. In Fig. [10] a sketch of the in- 
stallation is shown. Hot water cycled by a pump system 
is used to melt ice below the firn layer at 60 to 80 m 
depth to maintain a fresh water reservoir. An expanding 
cavity is formed with a diameter as large as 15 to 20 m. 
For IceCube and its predecessor AMANDA this tech- 
nique is used in connection with drilling at the South 
Pole since mid 1990s. If the well is used a second time 
a year later, a second cavern is formed at a deeper level. 

Having identified acoustic events arising from the 
2007/08 Rodriguez-Well, three other event clusters 
were found, two of them could be attributed to other 
IceCube Rodriguez-Wells from 2006/07 and 2004 to 
2006. The fourth event cluster turned out to be lo- 
cated at the probable position of the last AMANDA 
Rodriguez- Well used in the final two drilling seasons 
up to 2001. No documented coordinates could, how- 
ever, be found for that position. Available information 
about acoustic event clusters connected with Rodriguez- 
Wells is summarized in Table|2] As can be seen from the 
table and from Fig. |1 l(a)| the acoustic events from the 
two Rodriguez- Wells used only during one season are 
located at shallower depths than those from Rodriguez- 
Wells used twice. This is in agreement with expecta- 
tions from the sketch in Fig. [10] The former were seen 
to emit acoustic signals from regions of decreasing vol- 
ume around the well core and finally stopped, the older 
one in October 2008, the younger one in May 2009. As 
an example, Fig. [12] shows the time profile of refreez- 
ing for Rodriguez-Well 2007/08, which was used only 
during one season. In contrast to that, acoustic events 
are observed until today from the six and ten years old 
deeper wells (see Fig. 1 1 1(b)) . The mechanism of sound 
production in and around the Rodriguez-Well caverns is 
still under debate in particular for the older wells. 

4.3. Acoustic event simulation 

A simple acoustic transient event simulation is done 
by calculating the signal propagation times for the dis- 
tance d„ = tJ(x„ - x) 1 + (y„ - y) 2 + {z„ - z) 2 between 
source (e.g. IceCube hole at (x,y,z)) and sensors n = 
!>•••> imax with A?„ = d„/v s . The signal is transmitted 
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Figure 10: Section of the Camp Century well after a first and second 
season of operation (Schmidt and Rodriguez 1962) from [ 10]. 
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Figure 12: Re-freezing of Rodriguez- Well 2007/08 in the x-t plane. 
Acoustic signals are emitted from an decreasing area around the well 
core until they disappear in May 2009. The time is given in days 
relative to Jan. 1,2009. 
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Table 2: Positions of acoustic events from different Rodriguez-Wells. The errors are given by the rms values of the event distribution and the 
asterisk in the last column marks the date when our analysis stopped. 
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Figure 1 1 : Acoustic events distribution in (a) depth versus v-coordinate and (b) in v-coordinate versus time. Lines: fitted positions of AMANDA- 
RW (solid red), 04/05-05/06 IceCube-RW (dashed green), 06/07 IceCube-RW (dotted blue), 07/08 IceCube-RW (dashed-dotted magenta). The 
time is given in days relative to Jan. 1, 2009. 
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Figure 13: Difference between reconstructed event position and true 
event position for simulated events. Deviations seen in the upper re- 
gion (z > — 170 m) are caused by the depth dependent sound speed. 



from a random point inside a certain cylindrical volume 
(radius 2 m, depth 2000 m) around the source. 

Although knowing that the true IceCube hole diame- 
ter is about 0.7 m, we take into account the possibility 
that tension cracks might appear outside the hole bound- 
ing surface, suggesting a larger simulation radius. The 
reconstruction of events simulated with constant speed 
of sound and without considering attenuation effects im- 
plies an exact source localization, which is in contra- 
diction to the real data vertex results, where a specific 
spread of vertices around the source (Fig j7(b)) and a 
lack of events below and above a certain depth is vis- 
ible (Fig. |8(b)[ ). The major reason for misreconstruc- 
tion of events at shallow depth (-200 m < z < -1 m) 
is the depth dependence of the sound speed [2], which 
is therefore included in the simulation. Above 174.8 m 
and below 1 m depth, the parameterization 



-(262.379+199.833 



— I 5 — 1213.08 1— | S )ms _1 (3) 
ml I m I 



is used. Below a constant sound speed value of v s = 
3878 ms~' is assumed. Parameterization of Eq. [3]is 
obtained using a fit of in-ice sound speed data points 
These conditions are taken into account by inte- 
grating along the path between source and sensor using 
the parameterized value of v s at every point. Due to 
the absolute positions of the source and the sensors and 



their relative locations to each other, refraction effects 
are negligible as shown e.g. in Fig. 6 of reference 
and are thus neglected in this paper. Further improve- 
ments are achieved by using additional information on 
the acoustic pressure wave attenuation in ice. We apply 
the formula 



an 



Sq - do 



dn-dp 



(4) 



with the initial amplitude So at the distance do chosen 
to fit the real data and an attenuation length A = 300 m 
as measured for South Pole ice with SPATS JH. S^ is 
the corresponding signal amplitude at the sensor n. If 
the signal strength at a sensor is above ~ 300 mV (5.2cr 
above the noise level), a hit is triggered as in real data. 

Good agreement between reconstructed real and re- 
constructed simulated events is obtained below 170 m 
depth, where the localisation precision in the z coordi- 
nate is 25 cm (see FigfT3l. As expected, we observe 
a large influence of the depth dependence of the sound 
speed on reconstructions in the upper region of SPATS 
(between and 170 m depth), as one can also see in Fig. 
fT3~l The significant deviation of the sound speed from 
the constant value used in the reconstruction explains 
the spread of vertices seen in the real data (see Fig. 
|7(b)| i, whereas the direction of this smearing is caused 
by the detector geometry. Due to the strong attenuation, 
it is more difficult to observe deep events, which is well 
reproduced by the simulation. 

5. Estimated neutrino flux limit 

Although SPATS has not been built to measure a rel- 
evant neutrino flux limit, it is interesting to find out how 
sensitive a corresponding measurement could be, using 
data from this setup. In order to determine the number 
of events not connected to IceCube construction activi- 
ties in the sensitive region of SPATS, we omit the area 
of IceCube strings and the data from the drill periods. 
The area taken into account is indicated by the hatched 
area of Fig. |14(a)| Furthermore, we look at depths be- 
tween 200 and 1000 m, in the region of constant speed 
of sound, to avoid the smearing effect in the reconstruc- 
tion of acoustic event locations described in Section l4~3l 
In the 245 days of transient data taking (quiet period 2) 
we found no events in the defined region. This obser- 
vation is used to calculate an upper limit on the cosmo- 
genic neutrino flux. 

The effective target volume Fig. |14(b)| was calculated 
following the approach used for the first acoustic neu- 
trino limit estimate [11]. The neutrinos were assumed 
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to be down-going and to be uniformly distributed on a 
2n half sphere. The total cross-section is taken from 
a function derived by extrapolation of measured cross 
sections to higher energies II20I1 . which is valid for neu- 
trino energies above 10 5 GeV. Together with the interac- 
tion vertex, the direction (9, (f>) defines the plane of the 
acoustic pressure wave perpendicular to this direction. 
The sensor observation angle was then calculated rel- 
ative to this plane for each vertex. A number of 10 7 
events were simulated for neutrino energies E v from 
10 I8 eV to 10 22 eV. The energy £had of the hadronic 
cascade was assumed to be a constant fraction y = 0.2 
of the neutrino energy, i.e. £had = 0.2 E v . The contribu- 
tion of the cascade originating from the final state elec- 
tron in the electron-neutrino charged current reaction is 
omitted in the present model calculations, because its 
acoustic signal is expected to be small due to the LPM- 
effect [21]. The acoustic pressure P max was calculated 
with respect to observation angle and distance. We use 
the Askaryan lfl2ll model to calculate the acoustic sig- 
nal strength assuming a cylindrical energy deposition in 
the medium of length L and diameter d. No attempt was 
made to model angular sensitivity or frequency response 
of the sensor. A minimum threshold of ~ 300 mV, as 
in the real SPATS measurement, was applied. Using 
our estimate for the average SPATS sensor sensitivity 
(Sec. l3.3l >. this transforms to a necessary minimum pres- 
sure of ~ 70 mPa. At least five hits distributed over all 
four strings were required for an event to trigger. 

The procedure has a number of rather large associated 
uncertainties that are discussed in the following list: 

• The biggest uncertainty is related to the different 
predictions of the various thermo-acoustic models 
for the acoustic signal strength (see e.g. §Y~h . lfl8ln . 
We assume a 100 % uncertainty on this quantity. 
No absolutely calibrated measurement exists so far 
in any medium to fix that problem. 

• The Landau-Pomeranchuk-Migdal effect [ 19] adds 
an additional uncertainty by changing the cross 
sections of bremsstrahlung and pair-production at 
ultra-high energies. The effect elongates pref- 
erentially primary electron cascades from neu- 
trino interactions and diminishes expected acous- 
tic signals [21]. It starts to become important for 
hadronic cascades above 10 18 eV. Above 10 20 eV 
its influence is diminished by photo-nuclear and 
electro-nuclear interactions (see I122ll for a recent 
detailed discussion). 

• The uncertainty of about 40% in the absolute noise 
level determination leads to a corresponding uncer- 



tain trigger threshold. This makes the lower energy 
threshold for contributing neutrino interactions un- 
certain. 

• The angular efficiency loss of the single sensor 
channels is elaborated by use of a data set where 
two channels per sensor are available. With the 
threshold taken (70 mPa), an efficiency of > 99% 
is found for 99% of the azimuthal angular range. 
This leads to the conclusion that the azimuthal ef- 
ficiency loss is negligible in comparison to other 
error sources taken into account. 

• The 30% error on the sound attenuation length is 
of minor importance in comparison to the effects 
discussed above. 

These problems influence all acoustic (and partly ra- 
dio) neutrino limits given so far (see Fig. |14(c)} . In the 
present paper the acoustic signal is calculated with a 
model delivering signal values in the middle between 
extreme predictions. 

The observation of zero events inside the effective 
volume of SPATS gives an upper limit of A^bs = 2.44 
events at a Poissonian 90% confidence level l23"ll . The 
flux limit for an assumed trigger threshold of the mea- 
surement of 70 mPa is shown in Fig. |14(c)| as dashed 
curve together with cosmogenic neutrino flux predic- 
tions and results from other experiments. The grey and 
dark grey bands around the given limit indicate the un- 
certainties of this quantity due to the uncertainties dis- 
cussed above. The upper border of the grey shaded area 
can therefore be considered as a conservative neutrino 
flux limit derived from the SPATS data, provided that 
the assumptions made in this work for estimating the 
detector sensitivity hold. 

6. Summary and Outlook 

We presented an analysis of acoustic noise data 
recorded with the South Pole Acoustic Test Setup 
(SPATS) in the deep Antarctic ice at the geographic 
South Pole. We found the absolute noise level to be 
extremely stable over time. Its estimated magnitude 
of 14mPa in the deep ice, below 200 m, is compara- 
ble to the noise in the deep sea when weather condi- 
tions are calm Studies of transient noise in 
the SPATS data revealed the refreezing IceCube holes 
and Rodriguez-Wells as sources. The high quality of 
the data allowed us to study the refreezing processes as 
function of time in great detail. No transient acoustic 
signals in the deep ice were observed outside the instru- 
mented volume of IceCube at depths below 200 m. This 
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Figure 14: (a) Shown are the IceCube center (large black dot) with the construction area (black circle) and the SPATS center (large grey dot) with 
the sensitive area for acoustic detection (grey circle). In addition, strings with acoustic sensors (medium size dots) and all IceCube strings (small 
dots) are drawn in. To avoid events from IceCube construction, only the shaded area of the SPATS sensitive region is used in order to measure the 
neutrino flux limit, (b) Corresponding effective volume and (c) the neutrino flux limit of the 2009 SPATS configuration (70 mPa threshold, > 5 
hits per event). The dark gray band (50 to 100 mPa threshold) around the effective volume and around the limit considers uncertainties in absolute 
noise. The even broader light gray band includes additional uncertainties due to the choice of different acoustic models. Experimental limits on 
the flux of ultra high energy neutrinos are from ANITA II (Q, FORTE 0, GLUE (Hi, SAUND II [H, ACoRNE Q<|. For model reference see 

mm. 
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enabled us to derive a first upper limit on the flux of 
ultra-high energy neutrinos with an acoustic detector in 
glacial ice. 

SPATS is continuing to take data. An upgrade of the 
DAQ software to read out all sensor channels simulta- 
neously and to form a multiplicity trigger online, will 
increase the detector sensitivity. 
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